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We study the collective modes in the magnetic-field-induced spin-density-wave (FISDW) 
phases experimentally observed in organic conductors of the Bechgaard salts family. In 
phases that exhibit a sign reversal of the quantum Hall effect (Ribault anomaly), the coex- 
istence of two spin-density waves gives rise to additional long wavelength collective modes 
besides the Goldstone modes due to spontaneous translation and rotation symmetry break- 
ing. These modes strongly affect the charge and spin response functions. We discuss some 
experimental consequences for the Bechgaard salts. 

Introduction. - The organic conductors of the Bechgaard salts family (TMTSF)2X (where 
TMTSF stands for tetramethyltetraselenafulvalene) have remarkable properties in a mag- 
netic field. In three members of this family (X=C104, PFg, Re04), a moderate magnetic 
field of a few Tesla destroys the metallic phase and induces a series of SDW phases separated 
by first-order phase transitions [|^, 

According to the so-called quantized nesting model (QNM) P], the formation of the 
FlSDWs results from the strong anisotropy of these organic materials, which can be viewed 
as weakly coupled chain systems (the typical ratio of the electron transfer integrals in the 
three crystal directions is ta ■ tb : tc = 3000 : 300 : 10 K). The SDW opens a gap, but leaves 
closed pockets of electrons and/or holes in the vicinity of the Fermi surface. In presence of 
a magnetic field H, these pockets are quantized into Landau levels (more precisely Landau 
subbands). In each FISDW phase, the SDW wave vector is quantized, Qn = {^kp + NG, Qy) 
with integer, so that an integer number of Landau subbands are filled. (Here kp is the 
Fermi momentum along the chains, — e the electron charge, h the interchain spacing, and 
G = eHb/h.) As a result, the Fermi level lies in a gap between two Landau subbands, the 
SDW phase is stable, and the Hall conductivity is quantized: axy = —2Ne^/h per one layer 
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of the TMTSF molecules P, As the magnetic field increases, the value of the integer N 
changes, which leads to a cascade of FISDW transitions. 

A striking feature of the QHE in Bechgaard salts is the coexistence of both positive and 
negative Hall plateaus. While most plateaus are of the same sign, referred to as positive by 
convention, a negative Hall effect is also observed at certain pressures (the so-called Ribault 
anomaly) 0. 

We have recently explained the Ribault anomaly within the framework of the QNM by 
taking umklapp processes into account 0. Because of umklapp scattering, two linearly po- 
larized SDWs, with wave vectors Qtv = {^kp+NG, n/b) and Q-n = {^kp—NG, vr/6), coexist 
in the Ribault phase 0. The quantum Hall conductivity is quantized: axy = —2Ne'^/h. The 
integer N is negative (hence a negative Hall plateau) and corresponds to the SDW with the 
largest amplitude. Experimentally, = —2 and A^ = —4 are the most commonly observed 
'negative' phases. Our explanation of the Ribault anomaly differs from the one suggested 
by Zanchi and Montambaux ||^ by invoking the pressure dependence of umklapp scattering 
rather than the electron band structure. 

In this Letter, we study the long wavelength collective modes in the FISDW phases that 
exhibit the Ribault anomaly. The coexistence of two SDWs in these phases [0 gives rise to 
additional collective modes besides the Goldstone modes resulting from spontaneous rotation 
and translation symmetry breaking. We point out some analogies with phase modes in two- 
band or bilayer superconductors 0, [T^ and plasmon modes in semiconductor double-well 
structure [|lT|. While these modes are generally difficult to observe, collective modes have 



strong experimental consequences in SDW systems [|T2|. We discuss how the charge and spin 
response functions are affected in the Ribault phase. 

It should be pointed out that the long wavelength modes are not the only modes of interest 
in the FISDW phases. There also exist magneto- rotons at finite wave vectors (g^ = G, 2G, ...) 



T3| . Within the Zanchi-Montambaux scheme, Lederer has recently shown that these modes 



exhibit a different behavior in the Ribault phase []T1|. The effect of umklapp scattering on 
the magneto-rotons is not considered in this Letter. 

Mean-field theory. - In the vicinity of the Fermi energy, the electron dispersion law in the 
Bechgaard salts is approximated as 

E{kx, ky) = VF{\k^\ - kp) + tj_{kyb), (1) 

where k^ and ky are the electron momenta along and across the one-dimensional chains of 
TMTSF, and h = 1. In Eq. (|I]), the longitudinal electron dispersion is linearized in k^ 
in the vicinity of the two one-dimensional Fermi points ^kp, and vp = 2ata sm{kpa) is the 
corresponding Fermi velocity (a is the lattice spacing along the chains). The periodic function 
^±(w) = t±{u + 27r) describes the interchain hopping in a tight-binding approximation: 

tj_{kyb) = —2tbCos{kyb) — 2t2b cos{2kyb) — 2t3b cos{3kyb) — 2t4b cos{4:kyb) . (2) 

We neglect the electron dispersion in the third direction along the z axis, t^i, plays a crucial 
role in the Ribault phase since (together with the umklapp scattering strength) it determines 
the ratio of the amphtudes of the two SDWs. For t^f, = 0, both SDWs have the same 
amplitude (independently of the value of t^b)- A finite t^b lifts this degeneracy. This yields a 
negative QHE whenever sgnit^b) = sgn(t2b) [0- On the other hand, t^b does not affect directly 



the Ribault phase |15]. Its only effect is to change the critical value of umklapp scattering 



above which the Ribault phase becomes stable. As shown by Zanchi and Montambaux 
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a sufficiently strong t^b {tsb ~ 0.2t2fe) can stabilize a negative phase even in the absence 
of umklapp scattering. We cannot exclude that both umklapp scattering and dispersion 
law play a role in stabilizing the Ribault phase at low pressure. However, whatever the 
mechanism at work, umklapp scattering (not considered in the Zanchi-Montambaux theory) 
will always lead to the coexistence of two SDWs in the Ribault phase. 

The presence of a magnetic field along the z direction is taken into account via the Peierls 
substitution k = {k^, ky) —i'V — eA where A is the vector potential. The magnetic field 
also introduces the Zeeman coupling a^sH where a = +(— ) for up (down) spins, (/i^ is 
the Bohr magneton. We take the electron gyromagnetic factor g equal to two.) We consider 
forward and umklapp scatterings between electrons, with amplitudes g2 and 5^3 respectively. 
We do not consider backward scattering, since it does not play an important role in the 
QNM [§. 

In the FISDW phase, the electron spin density has a nonzero expectation value: 

(S'^(r)) = ^ m^ArCOs(0/3Ar)cOs(r-Q^7V + 6'/3Ar), 
/3=± 

{Sy{r)) = sin(0/3Ar) cos(r • Q/jat + 9^^), (3) 

/3=± 

where r = [x, y) is the spatial coordinate. Because of the Zeeman coupling with the magnetic 
field, the SDWs are polarized in the (x, ?/)-plane. 07v and 0_7v determine the direction of the 
spin magnetization, and 9^ and 9_n the position of the SDWs with respect to the crystal 
lattice. The ratio I7I = ui-m/^n of the amplitudes of the two SDWs in the Ribault phase 
is not precisely known. It depends on the detailed geometry of the Fermi surface (via ^4^) 
and increases with decreasing pressure. However, the stability of the Ribault phase against 

the formation of two helicoidal (i.e., circularly polarized) SDWs requires I7I ~ 0.5 0. 

The FISDW phases are characterized by the complex order parameters A^jv,Q defined by 

A.(r) = (V^L,,^(r)V^„,T(r)) = A;3^,.e--Q''-. (4) 

The operators il)^\{v) annihilate (create) electrons with spin a and momenta close to akp 
(a = ±). The order parameters are entirely determined by the mean value of the spin density 
(Eq. (^). In particular, |A^jv_|_| = |A^jv,-| = ±) for sinusoidal (i.e., linearly polarized) 
SDWs [|. 

First we consider the zero-temperature condensation energy /\E at the mean-field level 
using the quantum limit approximation (QLA), also known as the single gap approximation. 
This approximation is valid when vpG ^ T and consists in retaining only the gaps at the 
Fermi level, neglecting those opening above and below the Fermi level 0, p. Introducing 

A/3Ar,a = I/3N{,g2^l3N,a + g3^-/3N,~a), (5) 

we write the condensation energy as 



2En 



(6) 



Here i?o is an ultraviolet cutoff of the order of ta, and A(0) = Xj-nvpb the density of states 
per spin. The coefficients J„ = IniQyb = vr) are well known in the QNM. They depend on 
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the transverse dispersion law t±{kyb) and measure the degree of perfect nesting of the Fermi 
surface ^. 

Minimizing AE with respect to the order parameters Ap^^a, we find that the mean-field 
ground state corresponds to 6^ = —6-^ and (p^ = 4>-n- The latter relation shows that 
both SDWs have the same polarization axis. The condition 6^ = nieans that the two 

SDWs can be displaced in opposite directions without changing the energy of the system. 
This property is related to the pinning that would occur for a commensurate SDW. Indeed, 
for a single SDW with wave vector {2kF,n/b), the condition 9n + 6-n = becomes the 
usual pinning condition ^ = 0, where 9 is the phase of the SDW [O]. The degeneracy of 



the ground state results from rotational invariance around the z axis in spin space []T6| and 



translational invariance in real space. The latter holds in the FISDW phases, since the SDWs 



are incommensurate with respect to the crystal lattice [|T7 . 

According to the mean-field analysis, and in agreement with general symmetry consider- 
ations, we therefore expect two (gapless) Goldstone modes: a spin-wave mode corresponding 
to a uniform rotation around the z axis of the common polarization axis, and a sliding mode 
corresponding to a displacement of the two SDWs in opposite directions. 

Long wavelength collective modes. - Collective modes can be studied by expressing the 
partition function as a functional integral over bosonic fields describing spin fluctuations 



The standard mean-field theory is recovered from a saddle point approximation. Collective 
modes are obtained studying small (Gaussian) fluctuations around the saddle point |T^. The 



static mean-field order parameters Aa{r) (Eq. (^)) then become space- and time- dependent 
fluctuating variables Aa(r,r) = Aa{r) + 5Aa(r, r). Here r is an imaginary time. We use 
the Matsubara formalism and obtain real-time quantities by standard analytic continuation. 
Taking advantage of the decoupling of phase and amplitude modes in the long wavelength 
limit, we do not consider the latter. Order parameter fluctuations 6Aa{r, r) then correspond 
to fluctuations of the phase variables 9±n{t, t) and 0±7v(r, r). Skipping technical details pO| , 



we only quote the final result for the collective modes, restricting ourselves to longitudinal 
(i.e., parallel to the chains) fluctuations. 

We find two Goldstone modes with a linear dispersion law u = vpqx'- a sliding mode and 
a spin-wave mode corresponding to 

9N{qx,^^) = -9_N{q^,uj), 0Ar(g^,t^) = (p.^igx,^), (7) 

respectively. 9±n{Qx,^) and (f)±iy{qx,uj) are the Fourier transforms of the phase variables 
9±N{x,t) and (j)±N{x,t) (after analytic continuation to real time t). Eq. (|^) agrees with 
previous conclusions drawn from the mean-field analysis. Notice that the oscillations of the 
two SDWs are in-phase in the gapless spin-wave mode, and out-of-phase in the gapless sliding 
mode. 

We also find a gapped sliding mode, cj^ = v^ql + uf, and a gapped spin-wave mode, 
o;^ = Vpql + with: 



12 r 




3 + 572 


^2iV(0) 1 - r2 


InI-n 
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OJi = — ,^,„, — — z ; rrrr^ — , ^jJo 
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1-7^ 


^2A^(0) 1 - r2 


InI-n 


7^ 



where r = g^/g2 and 7 = A^n,+/An,+. 7 is related to 7 = A_n,+/An,+ (I7I = m_N/mN) 
by 7 = (7/Ar — tI^n) I {I^N — r'yiN)- In the Ribault phase, I7I ~ I7I [^. With no loss of 
generality, we choose Ap^^a and Ap^^a to be real. The gapped sliding mode corresponds to 
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a displacement of the two SDWs in the same direction (i.e., sgn[9N{x,t)] = sgn[9_]\f{x,t)]), 
while the gapped spin- wave mode corresponds to opposite rotations of the SDWs (i.e., 
sgn[(f)isf{x,t)] = — sgn[0_Ar(x, t)]). In particular, when I7I = I7I = 1 (a situation reached 
when I/tvI = |/-Ar|), gapped modes correspond to 9n = 9_n and 4>n = —^(p-N- Using the 
physical parameters of the Bechgaard salts, we find that ui and 002 are larger than the mean- 
field order parameters |A±jv±|, so that the gapped modes appear above the quasi-particle 
excitation gap (generally within the first Landau subband above the Fermi level). There- 
fore, we expect these modes to be strongly damped due to the coupling with quasi-particle 
excitations. 

There are some similarities between collective modes of the SDWs and phase modes 
occuring in two-band or bilayer superconductors [|, |T^. (73/(72 plays the same role as the 
ratio between the intraband (or intralayer) and interband (or interlayer) coupling constants. 
To some extent, there are also analogies with plasmon modes occuring in conducting bilayer 



systems [11|. While the corresponding phase modes in superconducting systems have not yet 
been observed, plasmon modes in semiconductor double-well structures have been observed 
recently via inelastic light scattering experiment |^ . 



Spectral functions. - Now we consider the spectral functions of the collective modes of the 
SDWs: the spin-spin correlation function Xyy (^5 ^) the optical conductivity cr(q, u) = 
{i / uj)W^^ {c[, uj) , where 11''*^* is the current-current correlation function. and x^yy ^''^^ the 
retarded parts of the imaginary time correlation functions 

n(r, r; r', r') = (jDw(r, r)jDw(r', r')), Xyyi^, r; r', r') = {Sy{r, T)Sy{r', r')). (9) 

jDW is the current along the chains carried by SDW fiuctuations. For real Af^^^a, the mean- 
field magnetization (Eq. (|)) is along the x axis so that Sy corresponds to transverse spin 
fluctuations. To lowest order in phase fiuctuations 



jDw(r,r) = — , 

nooT 1 + r7 

Sy{r,r) = 2^A^Ar,+ cos(r • Q/3Ar)0;37v(r,r). (10) 
P 

jdw is a function of 9±n, while Sy is a function of 4>±n- Thus, 11 is determined by the sliding 
modes, while Xyy is determined by the spin-wave modes. 

In the limit q = 0, the dissipative part of the conductivity is given by 



uj'/ .3(l-f) 4f 



Re[a{u)] = [ S{u)^—i^ + 6{u ± ^1)—%^ • (H) 



We have introduced the plasma frequency Up = ^Se^vp/b. Eq. (|TT|) satisfies the conductivity 
sum rule J^^duRe[a{u)] = uj^/A. Quasi-particle excitations above the mean-field gap do 



not contribute to the optical conductivity, a result well known in SDW systems [Q. Because 
both modes contribute to the conductivity, the low-energy (Goldstone) mode carries only a 
fraction of the total spectral weight. We obtain Dirac peaks at ±u!i because we have neglected 
the coupling of the gapped mode with quasi-particle excitations. Also, in a real system (with 
impurities), the Goldstone mode would broaden and appear at a finite frequency (below the 
quasi-particle excitation gap) due to pinning by impurities. In the clean limit, which is 
appropriate in (TMTSF)2X salts, the presence of impurities does not restore any significant 
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spectral weight to quasi-particle excitations above the mean-field gap [T^. Therefore, the 



fraction of spectral weight carried by the two modes is correctly given by Eq. (|TT]). By 
measuring the optical conductivity (j{ui), we can therefore obtain the ratio I7I ~ I7I of the 
amplitudes of the two SDWs. We have shown in Ref. [^] that I7I can vary between ~ 
and ~ 0.5 in the Ribault phase. Thus, the Goldstone mode can carry between ~ 100% 
and ~ 50% of the total spectral weight. When this fraction is reduced below ~ 50% (i.e., 

when I7I ~ 0.5), the Ribault phase becomes unstable against the formation of a helicoidal 
phase P]. The low-energy spectral weight is therefore also a measure of the stability of 
the sinusoidal Ribault phase against the formation of a helicoidal phase. For a helicoidal 
structure, we cannot distinguish between a uniform spin rotation and a global translation, 
so that there is only one type of modes. Thus, in the helicoidal phase, we find a gapless 
Goldstone mode and a gapped mode. The Goldstone mode carries no current, and all the 
spectral weight is pushed up above the quasi-particle excitation gap p]|. 



Due to the presence of SDWs, the spin-spin correlation function is not diagonal in 
momentum space, but has components X™ + aQ/3iv, q + a'Q/3'iv, uj) where a, a', /3, = ± 



(see Eqs. (PHIOD). q corresponds to the momentum of the spin- wave mode and tends to zero 
for long wavelength fiuctuations. We therefore consider the spectral function ImTrqX^^*, 
where Trq is a partial trace corresponding to a given spin-wave mode momentum q. In the 
limit g^, = 0, — 0, we obtain 



ImTrqX 



ret 



2n 


^N, + ^-N,+ \ 




\InI-n\ 




UJ2) 3(1 + r2 


+ — 

UJ2 


2(1 -r2) 




^2)2 l^^l 



(12) 



for a;,ga; > 0. C = I-n/In- Both spin-wave modes contribute to the spectral function. The 
spectral weight carried by the Goldstone mode diverges as 1/qx as expected for a quantum 
antiferromagnet |rP|. Eqs. (|TT]) and (|T^ predict that all the spectral weight is carried by 



the in-phase modes, i.e., the gapped sliding mode and the gapless spin- wave mode, whenever 
both SDWs have the same amplitude (I7I = I7I = |C| = 1). 

Conclusion. - We have shown that the Ribault anomaly in the FISDW phases of the Bech- 
gaard salts is characterized not only by a sign reversal of the QHE, but also by a rich 
structure of collective sliding and spin-wave modes. The presence of two SDWs in this phase 
gives rise to in-phase and out-of-phase collective oscillations in the long wavelength limit. 
The out-of-phase sliding and in-phase spin- wave modes are gapless (Goldstone modes). The 
in-phase sliding and out-of-phase spin-wave modes are gapped and are expected to appear 
above the quasi-particle excitation gap in the Ribault phase. Charge and spin response func- 
tions are strongly affected by the presence of these modes. In particular, we have shown that 
the low-energy (Goldstone) sliding mode carries only a fraction of the total spectral weight in 
the optical conductivity. By measuring the latter, we can obtain the ratio of the amplitudes 
of the two SDWs that coexist in the Ribault phase. The low-energy spectral weight is also 
a measure of the stability of the Ribault phase against the formation of helicoidal SDWs. 
Helicoidal SDWs have been predicted in Ref. and are expected to appear at low pressure 
(since I7I increases with decreasing pressure). We have pointed out some analogies with 
collective modes in two-band superconductors or bilayer systems |r^, |TT|. The common 



feature to all these collective modes is that they appear in a two-component system. To our 
knowledge, it is the first time that such modes are predicted in SDW systems. 
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Note that the long-range Coulomb interaction (not considered here) affects neither the 
transverse sliding modes sampled by optical spectroscopy nor the spin-wave modes. 
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